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We investigate the spin entanglement in the superconductor-quantum dot system proposed by 
Recher, Sukhorukov and Loss, coupling it to an electronic beam-splitter. The superconductor- 
quantum dot entangler and the beam-splitter are treated within a unified framework and the en- 
tanglement is detected via current correlations. The state emitted by the entangler is found to be 
a linear superposition of non-local spin-singlets at different energies, a spin-entangled two-particle 
wavepacket. Colliding the two electrons in the beam-splitter, the singlet spin-state gives rise to a 
bunching behavior, detectable via the current correlators. The amount of bunching depends on the 
relative positions of the single particle levels in the quantum dots and the scattering amplitudes 
of the beam-splitter. It is found that the bunching- dependent part of the current correlations is of 
the same magnitude as the part insensitive to bunching, making an experimental detection of the 
entanglement feasible. The spin entanglement is insensitive to orbital dephasing but suppressed by 
spin dephasing. A lower bound for the concurrence, conveniently expressed in terms of the Fano fac- 
tors, is derived. A detailed comparison between the current correlations of the non-local spin-singlet 
state and other states, possibly emitted by the entangler, is performed. This provides conditions for 
an unambiguous identification of the non-local singlet spin entanglement. 

PACS numbers: 03.67.Mn,73.50.Td,73.23.Hk 
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I. INTRODUCTION 



Ever since the concept of entanglement was 
introduced,* it has been at the heart of conceptual 
discussions in quantum mechanics The discussions 
have mainly concerned the non-local properties of entan- 
glement. Two entangled, spatially separated particles, 
an Einstein-Podolsky-Rosen 2 (EPR) pair, are correlated 
in a way which can not be described by a local, realistic 
theory, i.e. the correlations give rise to a violation of a 
Bell Inequality^ In optics, the non-local properties of 
entangled pairs of photons have been intensively inves- 
tigated over the last decadesi 6 i 7 i 8 Recently, the interest 
has turned to possible applications making use of the 
properties of entangled particles. Entanglement plays 
an important role in many quantum computation and 
information schemes^ with quantum cryptography 10 
and quantum teleportationii^ as prominent examples. 

Compared to optics, the investigation of entanglement 
in solid state systems is only in it's infancy. However, a 
controlled creation, manipulation and detection of entan- 
glement is a pre-requisite for a large-scale implementa- 
tion of quantum computation and information schemes, 
making it of large interest to pursue the investigation 
of entanglement in solid state systems. Considerable 
experimental 13 and theoretical^ progress has already 
been made in the understanding of entangled qubits im- 
plemented with Josephson junctions. 

For the entanglement of individual electrons, recently 
several important steps towards an experimental real- 
ization in mesoscopic conductors were taken. A scheme 
for entanglement of orbital degrees of freedom was pro- 
posed in Ref. 15], allowing for control of the entan- 
glement with experimentally accessible electronic beam- 
splitters>iSiii Moreover, several proposalsi^iSiiS for de- 



tecting entanglement via a violation of a Bell Inequality, 
expressed in terms of zero-frequency noise correlators^ 
have been put forth. Very recently, following a proposal 
by Beenakker et al [U, several works have discussed 
the possibility of electron-hole and post-selected electron- 
electron entanglement In particular, entanglement 
in the electrical analog of the optical Hanbury Brown 
Twiss effect* 2 ^ was investigated in a mesoscopic conduc- 
tor in the quantum Hall regime, transporting electrons 
along single edge-states and using quantum point con- 
tacts as beam-splitters. 23 Moreover, a scheme for energy- 
time entanglement^ has been proposed^ The conse- 
quences of dephasing for orbital entanglement have been 
investigatedi£*2L2& as well. 

Earlier proposals for electronic entanglement have 
been based on creating and manipulating spin en- 
tanglement, in normai22i22*2I as well as in normal- 
superconductin g 32 : 33 ! 34 systems. Spins in semiconduc- 
tors have been shown^ to have dephasing times ap- 
proaching microseconds, making spins promising candi- 
dates for carriers of quantum information. However, a 
direct detection of spin entanglement in mesoscopic con- 
ductors is difficult. The natural quantity to measure is 
the electrical charge current. To investigate spin cur- 
rent, one thus in principle has to convert the spin cur- 
rent to charge current via e.g. spin-filters. Although 
efficient spin-filters^ have very recently been realized 
experimentally, 37 there are considerable remaining ex- 
perimental complications in manipulating and detecting 
individual spins on a mesoscopic scale. In particular, to 
detect the entanglement by a violation of a Bell Inequal- 
ity, one needsiS two spin filters with independent and 
locally controllable directions to mimic the polarizers in 
optical schemes 

An alternative idea to detect spin entanglement was 
proposed by Burkard, Loss and Sukhorukov |30j and also 
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discussed qualitatively by Oliver, Yamaguchi and Ya- 
mamoto [18( . They proposed to use the relation between 
the spin and orbital part of the wavefunction, imposed 
by the antisymmetry of the total wavefunction under ex- 
change of two particles. A state with an antisymmet- 
ric, singlet spin wavefunction has a symmetric orbital 
wavefunction and vice versa for the spin triplet. When 
colliding the electrons in a beam-splitter, spin singlets 
and triplets show a bunching and anti-bunching behav- 
ior respectively. These different bunching behaviors were 
found to be detectable via the electrical current corre- 
lations, i.e. the properties of the orbital wavefunction 
were used to deduce information about the spin state. 
This approach was later extended to all moments of the 
current^ Moreover, it was recently further elaborated 
in Ref. [39l |. taking spin dephasing and non- ideal beam- 
splitters into account. 

In comparison to detecting spin entanglement via a vi- 
olation of a Bell Inequality, the approach of Ref. [sol ] 
however has a fundamental limitation. The antisymmet- 
ric spin singlet is an entangled state, while symmetric, 
triplet spin states are not neccesarily entangled. Consid- 
ering e.g. the standard singlet-triplet basis, only one of 
the three triplets | T4) + I IT), I TT) an d | ||) is spin en- 
tangled. However, all spin-triplet states, having the same 
symmetrical orbital wavefunction, give rise to the same 
anti-bunching behavior in the current correlators^ As 
a consequence, in contrast to a Bell Inequality test, the 
approach of Ref. j^lj can not be employed to distinguish 
between entangled and non-entangled triplet states. To 
be able to distinguish between different triplet state, one 
would need to consider more involved schemes, imple- 
menting in addition e.g. single spin rotations^ 

Despite this fundamental limitation, the approach of 
Ref. ;30( is due to its comparable simplicity still of 
interest for entanglers emitting non-local spin-singlets. 
However, the investigations in Ref. [30| were carried 
out assuming a discrete spectrum of the electrons and 
a mono-energetic entangled state incident on the beam- 
splitter. While giving a qualitatively correct picture of 
the physics, it does not quantitatively describe the situ- 
ation in a conductor connected to electronic reservoirs, 
where the spectrum is continuous and the entangled elec- 
trons generally have a wave-packet nature, i.e. the wave- 
function is a linear superposition of entangled electrons at 
different energies. 15 Moreover, the wavefunction in Ref. 
|30j | was not derived considering a specific entangler, it 
was instead taken to be an incoming plane wave with 
unity amplitude. This makes the calculated current cor- 
relations inapplicable to most of the entanglers consid- 
ered theoreticallyr2i22i2£ which operate in the tunneling 
regime and emit entangled states with a low amplitude. 

In this paper, we revisit the approach of detection of 
spin-singlet entanglement presented in Ref. [30[. The 
abovementioned shortcomings are bypassed by treating 
the entangler and the beam-splitter within a unified theo- 
retical framework. As a source of non-local spin-singlets, 
the superconductor-quantum dot entangler (see Fig. ^) 




FIG. 1: Schematic picture of the system. A superconductor 
(S) is connected, via tunnel barriers to two quantum dots (1 
and 2) in the Coulomb blockade regime. The dots are further 
coupled, via a second pair of tunnel barriers, to normal leads 
which cross in a forward scattering single-mode beam-splitter. 
The beam-splitter is characterized by scattering amplitudes 
r,t,r' and t' . On the other side of the beam-splitter, the 
normal leads are connected to normal electron reservoirs A 
and B. 



investigated in detail by Recher, Sukhorukov and Loss in 
Ref. [32| , is considered. Using a formal scattering ap- 
proach, the wavefunction of the electrons emitted from 
the entangler is calculated. It is found to be a linear su- 
perposition of pairs of spin-entangled electrons at differ- 
ent energies, a two electron wavepacket, similar to what 
was found for the superconducting orbital entangler in 
Ref. ^3 . The amplitude at each energy depends on the 
position of the single particle levels in the dots. Both the 
process where the electrons tunnel through different dots, 
creating the desired non-local EPR-pair, as well as the 
unwanted process when both electrons tunnel through 
the same dot, are investigated. In both cases the spin 
wavefunction is a singlet, preserving the spin-state of the 
Cooper pair tunneling out of the superconductor, how- 
ever the orbital states are different. 

The electrons emitted by the entangler are then col- 
lided in a beam-splitter and detected in two electronic 
reservoirs. Due to the singlet spin state, electrons tun- 
neling through different dots show a bunching behavior 
when colliding in the beam-splitter. Both the auto and 
cross correlations between currents flowing into the nor- 
mal reservoirs (but not the average current) depend on 
the degree of bunching. We find that the bunching is pro- 
portional to the wavefunction overlap of the two colliding 
electrons. This overlap depends strongly on the position 
of the single-particle levels in the dot, being maximal 
for both levels aligned with the chemical potential of the 
superconductors. The part of the current correlators sen- 
sitive to bunching is of the same magnitude as the part 
insensitive to bunching, making an experimental detec- 
tion of the spin-singlet entanglement feasible. 

The current correlators are independent of scattering 
phases and thus insensitive to orbital dephasing. How- 
ever spin dephasing generally leads to a mixed spin state 
with a finite fraction of triplets. Since the spin triplets 
have a tendency to anti-bunch, the spin dephasing re- 
sults in a reduction of the overall bunching behavior and 
eventually, for strong spin-dephasing, to a cross-over to 
an anti-bunching behavior. A simple expression for the 
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concurrence, quantifying the entanglement in the pres- 
ence of spin dephasing, is derived in terms of the Fano 
factors. 

For electrons tunneling through the same dot, the 
wavefunction is a linear superposition of states for the 
pair tunneling through dots 1 and 2. Both the cross- and 
auto correlators contain a two-particle interference term, 
sensitive to the position of the single-particle levels in the 
dots, however in a different way than the bunching depen- 
dent term for tunneling through different dots. In par- 
ticular, the correlators depend on the scattering phases, 
providing a way to distinguish between the two tunnel- 
ing processes by modulating e.g. the the Aharonov-Bohm 
phased Moreover, the phase dependence makes the cor- 
relators sensitive to orbital dephasing, while the spin part 
of the wavefunction is insensitive to dephasing. 



II. THE SUPERCONDUCTOR-QUANTUM DOT 
ENTANGLER 

A schematic picture of the system is shown in Fig. ^ 
A superconducting (S) electrode is connected to quantum 
dots (1 and 2) via tunnel barriers. The dots are further 
contacted, via normal leads to a controllable single-mode 
electronic beam-splitter— characterized by the forward 
scattering amplitudes r,t,r' and t' . The arms going out 
from the beam-splitter are connected to normal electron 
reservoirs A and B. 
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FIG. 2: Energy diagram of the entangler-beam-splitter system 
in Fig. Q A bias eV is applied between the superconducting 
reservoir, with chemical potential fis = 0, and the normal 
reservoirs A and B, with the same chemical potential fiN ~ 
—eV. There is only one spin-degenerate level of each dot, 
with energy ei and £2 respectively, in the energy range — eV 
to eV. The level width 7 is determined by the coupling to the 
normal reservoirs. The bias eV is taken to be much smaller 
than the superconducting gap, eV <C A, but so large that the 
broadened levels are well within the bias window, eV — \£j\ 3> 
7, J =1,2. 

We first concentrate on a description of the entangler, 
the superconductor-quantum dot part of the structure in 
Fig. ^ investigated in great detail in Ref. 32] . The 
entangler was also recently examined within a density 
matrix approach^ " The role of the beam-splitter is dis- 
cussed further below, after a discussion of the quantum 
state emitted by the entangler. To simplify our presen- 
tation we carry over the notation from Ref. |32( when 



nothing else is stated. 

An energy diagram of the superconductor-quantum 
dot-normal lead part of the structure is shown in Fig. 
121 A negative bias — eV is applied to the normal reser- 
voirs while the superconductor is grounded. The chem- 
ical potential of the superconductor is taken as a refer- 
ence energy, [is = 0, giving the chemical potential of 
both normal reservoirs una = ^nb = = —eV. Each 
dot 1 and 2 contain a single, spin-degenerate level in the 
energy range — eV to eV, with energy e\ and £2 respec- 
tively. The level spacing in the dots is assumed to be 
much larger than the applied bias, so no other levels of 
the dots participate in the transport. The temperature is 
much lower than the applied bias (but much larger than 
the Kondo temperature). 

The tunnel barriers between the dots and the super- 
conductor are much stronger than the tunnel barriers be- 
tween the dots and the normal leads. As a consequence, 
the broadening 7 of the levels in the dots (taken the same 
for both dots) results entirely from the coupling to the 
normal leads. The voltage is applied such that the entire 
broadened resonances are well within the bias window, 
i.e. eV — \sj\ 3> 7 with j — 1, 2. The quantum dots are 
in the Coulomb blockade regime, i.e. it costs a charg- 
ing energy U to put two electrons on the same dot. The 
ground state contains an even number of electrons in the 
lower lying levels, i.e. anti-ferromagnetic filling of the 
dots. 
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FIG. 3: Tunneling processes transporting two electrons from 
the superconductor to the normal leads. Process I, where the 
two electrons tunnel through different dots, one through dot 
1 and one through dot 2, creates the wanted EPR-pair. Pro- 
cess II, where the two electrons tunnel through the same dot, 
either both through dot 1 or both through dot 2, is unwanted. 

The transport takes place as Cooper pairs tunnel from 
the superconductor, through the dots and out into the 
normal leads. Due to the dominating tunnel barrier at 
the dot-superconductor interface, one pair that tunneled 
onto the dots leaves the dots well before the next pair tun- 
nels. There are two distinct possibilities for the Cooper 
pair to tunnel from the superconductor to the normal 
leads, shown in Fig. [3J 

• I, the pair splits and one electron tunnels through 
each dot, 1 and 2. 

• II, both electrons tunnel through the same dot, 1 
or 2. 

It was shown in Ref. [3^J that under the conditions stated 
above, all other tunneling processes could be neglected. 
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The process I creates the wanted EPR-pair, a spin singlet 
state with the two electrons spatially separated. How- 
ever, in an experiment one can not a priori exclude the 
second, unwanted process, II. One thus has to investigate 
process II as well, to provide criteria for an unambiguous 
experimental identification of emission of EPR-pairs. 

The first process, I, with the two electrons tunnel- 
ing through different dots, is suppressed below the sin- 
gle particle tunneling probability squared, since the two 
electrons have to leave the superconductor from two 
spatially separated points, i.e. effectively breaking up 
the Cooper pair. The tunneling amplitude for a bal- 
listic, three dimensional superconducto»22i2i is Aq oc 
exp(— d/^)/(kFsd) where d is the distance between the 
superconductor-dot connection points, kps the Fermi 
wave number in the superconductor and £ the supercon- 
ducting coherence length. This amplitude is in general 
larger for lower dimensional^ and disordered^i^ super- 
conductors. An investigation of the dependence of Aq on 
the geometry of the contacts to the superconductor was 
performed in Refs. 001- We point out that ways to 
avoid the suppression due to pair breaking by means of 
additional dots have been discussed in a similar context 
in Ref. |4jJ. However, since more dots complicate the 
calculation as well as the experimental realization, we 
consider the simpler geometry in Fig. ^ 

The second process, II, with both electrons tunneling 
through the same dot, is suppressed by the Coulomb 
blockade in the dots, as 1/U. In addition, there is a 
process which avoids double occupancy of the dots but 
instead requires a pair breaking, leading to suppression 
of the order 1/A. Together, this gives an amplitude 
Bq cx (1/U + 1/Air). The exact expression for the con- 
stants Aq and Bq in terms of tunnel amplitudes between 
the dots and the superconductor and the dots and the 
leads can be found in Ref. 32], for our purposes these 
expressions are not necessary. 

We point out that possibles candidates for experimen- 
tal realization of the proposed system are the extensively 
investigated^ heterostructures with semiconductors con- 
tacted to metallic superconducting electrodes. Electron 
transport through double dots in semiconductor systems 
have been recently been reviewed^ with an emphasis on 
experimental advances. 



III. THE WAVEFUNCTION OF THE 
SPIN-ENTANGLED ELECTRONS. 

To calculate the wavefunction of the electrons emitted 
from the superconductor-quantum dot entangler, we em- 
ploy the formal scattering theory 48 with the Lippman- 
Schwinger equation expressed in terms of the transfer 
matrix (T- matrix). The total Hamiltonian of the system 
can be written as H = Hq + Ht, where Hq is the Hamil- 
tonian of the superconductor, the quantum dots, and the 
normal leads. The perturbation Ht describes tunnel- 
ing between the superconductor, dots, and leads. The 



exact many-particle state \^f) satisfies the Schrodinger 
equation (E — H)\^f) = 0. In the absence of a per- 
turbation, Ht = 0, the system is in the ground state 
|0) = |0)s|0)d|0)jv, with different chemical potentials, 
fis = 0, and [In = —e.V. The perturbation Ht causes 
the electrons to tunnel from the superconductor, via the 
quantum dots, to the normal leads. 

We use the local nature of the tunneling perturbation 
and take the formal scattering approach to the prob- 
lem. According to this approach the state can be 
obtained by solving the Lippman-Schwinger equation in 
Fock-space 

|*) = |o) + 6(0)^1*), (i) 

where the retarded operator G(E) = [E— Hq+iO] -1 gives 
a state describing particles going out from the scattering 
region. Note that the total energy of the ground state |0) 
is E = 0. The formal solution of Eq. QJ can be written 
as 

|*) = |0>+G(0)T(0)|0), (2) 

where 

oo 

T(E) = H T + H T ^[G(E)H T ] n (3) 

n=l 

is the T-matrix. One then inserts a complete set of many 
body states 1 = J2n W)(N\ with \N) the eigenbasis of 
the Hamiltonian Hq, i.e. the basis of Fock-states of elec- 
trons and quasiparticles in the leads, dots and supercon- 
ductor respectively. The quantum number N collectively 
denotes the energies, spins, lead and dot indices etc of the 
individual particles. The eigenenergy of the state \N), i.e. 
the total energy of the individual particles, is En. This 
gives an expression for the state 

l*> = 10) -Y,Y^M lN){NlT{m - (4) 

N N 

In the system under consideration, all relevant matrix 
elements^ (jV|T(0)|0) are analytic in the upper part of 
the complex energy plane. As a consequency, in the in- 
tegration over energies of the individual particles in \N), 
the pole arising from the denominator En — iO can be 
replaced by a 5 (En) -function, imposing a total energy 
En = 0, equal to the chemical potential energy of the 
superconductor. This gives the wavefunction 

|*) = |0) - 2 m Y / S(E N )\N)(N\T(0)\Q). (5) 

JV 

It was shown in Ref. [22| , that under the conditions stated 
above and to lowest order in coupling between the su- 
perconductor and the dots, the operator T creates from 
the vacuum |0) a two-electron spin-entangled state. As 
pointed out above, depending on the relation between the 
amplitudes Aq and Bq, the transport of the two electrons 
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through the same (process II) or different (process I) dots 
dominates. Below we consider for simplicity only the lim- 
iting cases, where either I or II is completely dominating, 
however our analysis can straightforwardly be extended 
to a situation where they are of comparable strength. 



A. Electrons tunneling through different dots. 

We first consider process I, when the amplitude for 
tunneling through different dots is much larger than the 
amplitude to tunnel through the same dot. This creates 
the desired EPR-pair, a non-local spin-entangled pair of 
electrons. The quantities in this limit are denoted with a 
I. The wavefunction for two spin-entangled electrons at 
energies E\ and E 2 is 

\E u E 2 )j = i= [6^(^)6^ (^ 2 )_ bl^E^E, 

(6) 

where the operator b\ a (E) creates an outgoing (from the 
dots towards the beam-splitter) electron plane wave with 
spin a =t, [ and momentum k{E) — kp + E/hvp in the 
normal lead I = 1,2. Here kp and vp is the Fermi wave 
number and velocity respectively, same for both normal 
leads. The amplitude for this process was found in Ref. 
321 to have a double-resonant form 



(E 1 ,E 2 \T(0)\0) I = 



iA Ql /(irV2) 



{Ei +ei- il/2){E 2 + e 2 - h/2) ' 

(7) 



With this we are able to obtain the asymptotics of the 
outgoing spin-entangled state. For doing so we substitute 
Eq. J7J) into Eq. © and find 



If/ 



eV 



eV 



dEA(E) [&I T (£)^ (-25) - 6 f u (£)4 T (-E)]\0) 



with 



A(E) 



{E + ei- ij/2){-E + e 2 - iy/2) 



(8) 



(9) 



i.e. A(E) = (-i-jrV2)(E,-E\T(0)\0)i. This state is the 
sum of the unperturbed groundstate and an entangled, 
two electron state. The entangled state is a linear super- 
position of spin singlets at different energies, an entangled 
two-particle wavepacket. The singlet spin-state results 
from the singlet state of the Cooper-pair, conserved in 
the tunneling from the superconductor. Moreover, the 
two electrons in each singlet have opposite energies E 
and —E (counted from fig — 0), a consequence of the 
Cooper-pairs having zero total energy with respect to 
the chemical potential of the superconductor. 

Several important observations can be made regarding 
the state in Eq. (JSJ). First, the properties, including the 



two-particle wavepacket structure, can be clearly seen by 
writing the wavefunction in first quantization. Introduc- 
ing |l,i£)(g)| j) for bWE)\0), the properly symmetrized 
wavefunction is given by (omitting the ground state |0)) 

reV 

I*/} = / dEA(E){\l,E) y ,\2,-E) v + %-E) ll \l,E) v ) 

J-eV 

<8> (I T) M | iu) - | 1%) (io) 

with fi, v the particle index. The coordinate dependent 
wavefunction ^i(x fi ,x l/ ) = x v \^i) can then be writ- 
ten (x = at the lead-dot connection points) 



*i(x», x v ) = Mx^xu) (A* A* + xlxl) (xlxi - xlxl) 



(11) 



with 



ip(xu,x u ) 



2Tri^Ao 



2e — ij 

x exp [ik F (x u + Xp) - i(e - ij/2)\x u - x^/Hvf] 

(12) 

where for simplicity the case with energies E\ — e 2 = e 
is considered. To arrive at Eq. Hill we hrst introduced 
the wavefunctions (xi\E) ~ exp[ik(E)xi], I = ji,v, the 



spin spinors (xi\ |) = xli{xi\ I) = Xi an( i l ne orbital 
spinors (xi\l) = A, 1 , (xi\2) = Xf and then performed the 
integral over energy. The orbital spinors describe the 
wavefunction in the space formed by the lead indices 1 
and 2, a pseudo-spin space, as discussed in Ref. 0. We 
note that the beam-splitter, discussed below, only act in 
the orbital 12-space (i.e. spin independent scattering). 
Moreover, it is the property of the state in 12-space that 
determines the current correlators discussed below. 

As is clear from Eq. (|llfl , the state is a direct product 
state between the spin and orbital part of the wavefunc- 
tion. The spin state is antisymmetric under exchange 
of the two electrons, a singlet, while the orbital state is 
symmetric, a triplet. The probability to jointly detect 
one electron at x^ in lead 1 and one at x v in lead 2 decay 



exponentially with the distance \x^ 



an effect of the 



two electrons being emitted at essentially the same time 
(separated by a small time h/A) to points x^ — and 
x v = respectively. Note that the state |*/), a station- 
ary scattering state, does not describe wave packets in 
the traditional sense with two electrons moving out from 
the dots as time passes (as a solution to the time inde- 
pendent many particle Schrodinger Equation, 1^/) has a 
trivial time dependence). To obtain such a wavefunction, 
one must break time translation invariance by introduc- 
ing a time dependent perturbation, e.g. a variation of 
the tunnel barrier strength or dot-level energies in time. 

In this context it is worth to mention that such a time 
dependent wavefunction was recently considered by Hu 
and Das Sarma 1491 for a double-dot turnstile entangler. 
However, in Ref. |49| . the entangled wavefunction was 
not derived from a microscopic calculation but merely 



6 



postulated. The wavefunction had an amplitude of or- 
der unity (no tunneling limit) and contained a double 
integral over energy. This is different from our wavefunc- 
tion in Eq. JSJ) and moreover gives rise to a qualitatively 
different result for the currents as well as the current 
correlators studied below. 

Second, the entangled state in Eq. © has just the 
same form as the pair-splitted state obtained in the 
normal-superconducting system of Ref. 15], where a 
scattering approach based on the Bogoliubov-de Gennes 
equation was used. This shows rigorously that the effect 
of the strong Coulomb blockade, prohibiting two elec- 
trons to tunnel through the same dot, can be incorpo- 
rated in a scattering formalism by putting the amplitude 
for Andreev reflection back into the same dot to zero. 
From this observation it follows that the rest of the calcu- 
lation in the paper where the state in Eq. (JSJ is employed 
could in principle be carried out strictly within the scat- 
tering approach^ to the Bogoliubov de Gennes equation. 
However, in such a calculation the entanglement is not 
directly visible, which makes the interpretation of the re- 
sult difficult. Instead, below we work directly with the 
state in Eq. © . 

Third, it is also interesting to note the close connection 
between the emission of a Cooper pair and the process 
of spontaneous, parametric down-conversion^i of pairs 
of photons investigated in optics, where a single photon 
from a pump-laser is split in a non-linear crystal into two 
photons. From the point of view of the theoretical ap- 
proach, expanding the outgoing state in a ground state 
and, to first order in perturbation, an emitted pair of 
particles, is similar to the work in e.g. Ref. [5^]. The re- 
sulting state, Eq. iJHJ, is a spin singlet, while a state with 
polarization entanglement is, under appropriate condi- 
tions, produced in the down-conversion process (type II). 
Moreover, the emission of the two electrons is "sponta- 
neous" , i.e. random and uncorrelated in time, just in 
the same way as for the down-converted photons. One 
can also point out the maybe less obvious relation that 
the two electrons emitted from the superconductor carry 
information about the phase of the superconducting con- 
densate, just as the two photons carry information of the 
phase of the field of the pump-laser. A coherent superpo- 
sition of states of pairs of electrons emitted from different 
points of the superconductor, can give rise to observables 
sensitive to the difference in superconducting phase be- 
tween the two emission points, as was demonstrated in 
Ref. (lEf, This has its analog in the photonic experi- 
ment with a single, coherent laser pumping two separate 
non-linear crystals, presented in Ref. |53| . 



B. Electrons tunneling through the same dot. 



E% in lead j is 



\Ei,E 2 ) u = -Llb^E^Ei) - fct^Jftt^p), 

(13) 

The amplitude for this process, (E\,E2\T(Q)\0)u, was 
found in Ref. [32j to have a single resonant form, different 
from Eq. Q, 



(E U E 2 \T{Q)\Q) H 

1 



iB 
n2V2 



E! 



ij/2 



i 7 /2 



(14) 



Here, for simplicity the two dot-superconductor contacts 
are taken to be identical. Since the superconductor is a 
macroscopically coherent object, the total state is a linear 
combination of the states corresponding to two electrons 
tunneling through dot 1 and dot 2. To obtain the asymp- 
totics of the outgoing spin-entangled state, we substitute 
Eq. {nj into Eq. (J5J and find 



eV 



dE 



eV 



B 2 (E)bl 1 (E)bl l (-E) 



with 



B ( £j - in 12) 



(25- 



*7/2)(-25- 



i 7 /2) 



(15) 



(16) 



i.e. Bj(E) — (—iir2y/2)(E,—E\T(0)\0}n. Arriving at 
Eq. (JTSJ we used the property B(-E) = B(E) and the 
anti-commutation relations of the fcrmionic operators. 

This state is a linear superposition of the states for two 
electrons tunneling through the same dot. Comparing to 
the state i) in Eq. (jllfl for the two electrons tunnel- 
ing through different dots, we can make the following 
comments: (i) Just as \^>i), the wavefunction \^u) in 
first quantization is a product of a an orbital and a spin 
wavefunction. The spin wavefunction is, as for |\E'j), a 
singlet xjiX'l ~ XaXl- The orbital wavefunction for the 
simplest situation si = £2 is however proportional to 
A* A* + A^A^, one of the Bell states, an orbitally entan- 
gled state, (ii) The state \^fn) is the same that would be 
obtained within scattering theory (as was shown in Ref. 
0]), taking Bj(E) to be the effective Andreev reflection 
amplitude at dot j and assuming no crossed Andreev re- 
flection between the dots, i.e. zero probability for an 
incident electron in lead 1 to be back-reflected as a hole 
in lead 2 and vice versa. 

With the state in Eq. (|15f) and the state for two elec- 
trons tunneling through different dots, in Eq. JHJ), we are 
in a position to analyze the transport properties. 



We then turn to process II, when the amplitude for 
tunneling through the same dot is much larger than the 
amplitude to tunnel through different dots. The wave- 
function for two electrons to tunnel to energies E\ and 



IV. CURRENT CORRELATORS 

The two electrons emitted from the dot- 
superconductor entangler propagate in the leads 1 
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and 2 towards the normal reservoirs A and B. As shown 
in Fig. 01 the two normal leads are crossed in a single 
mode reflectionless beam-splitter. The beam-splitter is 
characterized by a spin- and energy independent unitary 
scattering matrix connecting outgoing and ingoing 
operators as 



r a 

t r' 



(17) 



where the subscript a — A, B denotes towards what 
reservoir the electron is propagating. The electrons are 
then detected in the normal reservoirs A and B. 

We point out that beam-splitters without backscat- 
tcring are not easily produced experimentally, thus in a 
more detailed model one should also take the effect of 
back scattering into account. Several aspects of back 
scattering were recently investigated b y B urkard and 
Loss extending the model in Ref. [3(|. Although 
back-scattering can be incorporated in our model as well, 
this would complicate the calculations and make the re- 
sult less transparent. Below, we instead neglect the effect 
of backscattering (however, some qualitative aspects are 
discussed below, in the context of the Aharonov-Bohm 
effect), pointing out that result below holds rigorously 
only in the case with negligible back-scattering at the 
beam-splitter. 

We also note that in the typical system of interest, 
with a lateral size L in the micrometer range, the energy- 
dependent part of the phase ~ Lj/Hvf picked up by the 
electrons when propagating in the leads is negligiably 
small. Energy independent phases due to propagation 
can be included in the scattering amplitudes of the beam- 
splitter. 

The properties of the electrons emitted by the en- 
tangler are investigated via the current and the zero- 
frequency current correlators. The electrical current op- 
erator in lead a is given by*^ 



e 
h 

x J2l b aAE)b a «(E')-aUE)a a<7 (E')], (18) 

a 

where a' a(J (E) creates an electron plane wave incoming 
from the normal reservoir a with spin a =f, J, and mo- 
mentum k(E). The averaged current is given by 

I a = (4), (19) 

where (...) = The zero- frequency correlations 

between the currents in the leads a and f3 are 

S aP = ( ' dt(Ai a (t)Aip(0) + Alp(0)At a (t)), (20) 



where AI a (t) = I a (t) — I a is the fluctuating part of the 
current in lead a. We study the two cases with elec- 
trons tunneling through different dots and the same dot 
separately. 



V. TUNNELING THROUGH DIFFERENT 
DOTS. 

For electrons tunneling through different dots, the 
question is how the degree of spin-singlet entanglement 
is reflected in the current and current correlators. The 
averaged current, evaluated with the state in Eq. 
©, becomes 



li = T J dE\A(E)\ 



(21) 



same for both a = A and B. Since the two resonances 
£i and £2 are well within the voltage range, i.e. eV — 
|ei|, eV — \e 2 \ ^> 7, we get the current 



li = 



2e 



4^o 12 
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h (e 1 +e 2 ) 2 + 7 2 ' 



(22) 



just the same expression as in Ref. [32|, where the leads 
of the entangler were contacted directly to the normal 
reservoirs (no beam-splitter). The current is maximal 
for an asymmetric setting of the resonances e± = — £2- 
This two-particle resonance reflects the fact that the two 
electrons in the Cooper pairs are emitted at opposite en- 
ergies with respect to the superconducting chemical po- 
tential. The current contains no information about the 
entanglement of the emitted state. In fact, the same cur- 
rent would be obtained by considering a product state of 
one electron in lead 1 and one in lead 2, independent of 
their spins. 

To obtain information about the entanglement, we 
turn to the current correlators. Inserting the ex pres sion 
for the state into Eq. following Ref. [54|, we 

get the expressions for the auto-correlations 

S'aa = S'bb = ^ dE{[l + 2RT]\A{E)\ 2 

+2RTA(E)A*(-E)} (23) 
as well as the cross-correlations 

S I ab = S i ba = ^ J*\e{[T 2 + R 2 ]\A(E)\ 2 

-2RTA(E)A*(-E)} , (24) 



where R = 



and T = \t\ 2 = \t'\ 2 = 1 - R. We 



note that the total noise S 1 of the current flowing out of 
the superconductor is twice the Poissonian, i.e. 

S 1 = S{ B + S BA + S AA + S BB = Ae(l{ + /£), (25) 

describing an uncorrelated emission of pairs of electrons. 
This result, an effect of the tunneling limit, is different 
from the one in Ref. [3(J where an entangled state with 
unity amplitude was considered and the total noise was 
found to be zero. 

It is clear from the calculation that the second term 
in Eqs. I|23|l and (|24|l depends directly on the symmetry 
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properties of the orbital wavefunction, and thus, due to 
the anti-symmetry of the total wavefunction, indirectly 
on the symmetry properties of the spin wavefunction. For 
a spin-triplet state the last term in Eqs. 1231) and 
(|24|l would have opposite sign. Since all the three possi- 
ble triplets, with spin wavefunctions xjiXl + XuXvi XaXl 
and xjiXl have the same anti-symmetric orbital wave- 
function (A* A 2 — A 2 A* for e\ = £2) they give rise to the 
same noise correlators. As a consequence, performing a 
noise correlation measurement, one can only distinghuish 
between spin-singlets and spin-triplets, but not between 
entangled x\xl + xjtXl and non-entangled xlxl, xltXi 
spin-triplets. This was pointed out already in Ref. [20] . 
We note that it is possible to distinguish between the dif- 
ferent triplets in a more advanced beamsplitter scheme, 
using controlled single spin rotations via a e.g. local 
Rashba interaction^ Such a scheme is straightforwardly 
included into our theoretical treatment, however, it de- 
mands a more involved experimental setup and is there- 
fore not considered here, we restrict our investigation to 
the simplest possible system. 

To investigate the properties of the current correlators 
in detail, the remaining integral over energy in Eqs. (IL'MI) 
and lj24")l is carried out, giving 



J dEA{E)A*{-E) = — 



47t|A | 2 7 3 



[(ei-e 2 ) 2 +7 2 ][(ei+£2) 2 +7T 

(26) 

This shows that, unlike the current, the noise is sensi- 
tive to both the difference and the sum of the dot energy 
levels. We note that the integral of A(E)A*(— E) is man- 
ifestly positive and smaller than the integral of ^(E 1 )! 2 
for all £i, £2 except for e\ = £2, when they are equal. 

From these observations we can draw several conclu- 
sions and compare our results to the results in Ref. [30l| : 

(i) The second term in Eqs. I|23|) and (|24|l . dependent 
on the orbital symmetry of the wavefunction, leads to a 
suppression of the cross-correlation but an enhancement 
of the auto-correlation. This is an effect of the bunching 
behavior of the spin-singlet, i.e. the two electrons show 
an increased probability to end up in the same normal 
reservoir^ For a symmetric beam-splitter, R = T = 
1/2 and aligned dot-levels £1 = £2, the cross-correlations 
are zero (to the leading order in tunneling probability 
considered here) . This is a signature of perfect bunching 
of the two electrons. 

(ii) The last term in Eqs. I)23|) and l|24|) is proportional 
the spectral overlap J dEA(E)A*(— E). The spectral 
overlap physically corresponds to the overlap between the 
wavefunctions of the two electrons colliding in the beam- 
splitter. For single-particle levels at different energies 
£1 7^ £2, the spectral amplitudes of the emitted electrons 
are centered at different energies and consequently^ the 
Pauli principle responsible for the bunching is less effi- 
cient. 

It is important to note that the last term in Eqs. ll'-il) 
and Q24)). dependent on the bunching, generally is of the 
same magnitude as the first term. We emphasize that 



this result is qualitatively different from what was found 
in Ref. [3(| , where the bunching dependent part of the 
current correlator was proportional to a Kronecker delta- 
function in energy, a consequence of considering a dis- 
crete spectrum. Our result clearly shows that it should 
be experimentally feasible to detect the bunching, and 
thus demonstrate that spin singlets are emitted from the 
entangler. We note that the same qualitative result was 
found in Ref. |49|. 

(Hi) The cross correlations are positive for any trans- 
parency of the beam-splitters (note that R 2 +T 2 > 2RT). 
This is different from the result in Ref. j2jj, where 
negative cross correlations were predicted. The nega- 
tive correlations are again a result of the unity ampli- 
tude of the incoming entangled state considered in Ref. 
[30|. In this context, we point out that positive cross- 
correlations have been predicted in several few mode^ 
and many mode^ normal-superconductors hybrid sys- 
tems as well as purely normal systems in the presence 
of interactions S In several of these cases, the positive 
correlations were explained with semiclassical models. 
Thus, the presence of positive correlations can not itself 
be taken as a sign of spin-entanglement. 

We point out that the expression for the energy de- 
pendent integrand of the cross correlators in Eq. I|24|) 
can be understood in an intuitive way, by considering 
the elementary scattering processes contributing to the 
noise, shown in Fig. Let us consider the probability 




FIG. 4: Elementary scattering processes (shown at the beam- 
splitter) contributing to the cross correlators S AB . The two 
processes (a) and (b) transport a pair of electrons |, E) and 
I l,—E) from the superconductor to the reservoirs A and B 
respectively. The two processes, having the same initial and 
final state, are indistinguishable and their amplitudes must be 
added. The correlator S ! AB is proportional to the integral over 
energy of the (energy dependent) joint detection probability. 

for the two electrons emitted from the superconductor to 
end up, one with spin up and energy E in reservoir A 
and the other with spin down and energy —E in reser- 
voir B. There are two paths the electrons can take from 
the superconductor to the reservoirs: (a) the electron 
with spin up and energy E via dot 1 and the electron 
with spin down and energy —E via dot 2. This pro- 
cess has an amplitude tt'A(E) (b) the electron with spin 
down and energy —E via dot 1 and the electron with 
spin down and energy —E via dot 2. This process has 
an amplitude rr'A(-E). Since the two processes have 
the same initial and final states, they are indistinguish- 
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able and their amplitudes must be added. This gives 
together the energy dependent joint detection probabil- 
ity ~ \tt'A{E)+rr'A(-E)\ 2 = T 2 \A(E)\+ R 2 \A{~E)\ 2 + 
rr l t*t*A(E)A*{-E)+r*r'*tt'A{-E)A*{E). In analogy 
to the noise correlators for the entangler with energy in- 
dependent tunneling probabilities in Ref. .15], it is found 
that the noise correlator S AB is simply proportional to 
integral over energy of the joint detection probability. 
Using that the integral in Eq. i|24|) goes from — eV to 
eV and that the unitarity of the scattering matrix in Eq. 
(|17|l gives rt* + t'r * = 0, we get the expression in the 
integrand in Eq. 

For the auto-correlation, a similar interpretation in 
terms of probabilities for two-particle scattering pro- 
cesses only is not possible, one also has to consider sin- 
gle particle probabilities. Formally, this is the case since 
auto-correlations contain exchange effects between the 
two particles scattering to the same reservoir. 




(e r e 2 )/y (e r e 2 )/y 



FIG. 5: The Fano factor for the cross-correlations Fab = Fba 
(left panel) and auto-correlations Faa = Fbb (right panel) 
as a function of the normalized energy difference (ei — £2)/"/ 
for various beam-splitter transparencies. 



A. Fano factors 

A quantitative analysis of the current correlators is 
most naturally performed via the Fano factors F a p = 
S a p/(2ey/I a Ip)- The Fano factor isolates the de- 
pendence of the noise on various parameters, not al- 
ready present in the current. For the cross- and auto- 
correlations respectively, we have 

Fab = Fh A = T 2 + R 2 - 2RT\H{e l - e 2 )\ 2 (27) 

and 

Faa = F i bb = 1 + 2RT + 2RT\H(e 1 - e 2 )\ 2 (28) 
where 



We note that only the last terms in Eqs. l(2Tjl and f2"%|l 

depend on the energies £1 and £2 of the levels in the 
dots. The Fano factor as a function of energy difference 
£1 — £2 is plotted for several values of transparency of 
the beam splitter in Fig. 03 For the cross-correlators, 
the Fano factor has a minimum for the two resonant lev- 
els aligned, S\ = £2. The value at this minimum de- 
creases monotonically from 1 to when increasing the 
transparency T of the beam-splitters from to 0.5 (the 
Fano factor for transmission probability T is the same 
as for 1 — T). Thus, for a completely symmetric beam- 
splitter, T = R = 0.5, the Fano factor is zero. This 
corresponds to the case of perfect bunching. For the 
auto-correlators, the picture is the opposite. The Fano 
factor has a maximum for the two resonances aligned, 
£1 = £2. The value at this maximum increases monoton- 
ically from 1 to 2 when increasing the transparency T of 
the beam-splitters from to 0.5. Thus, for a symmetric 
beam-splitter, T = R = 0.5, the Fano factor is now two. 



B. Decoherence 

Considering the robustness of the bunching behavior, 
an important observation is that the Fano factors in Eqs. 
(|28|l and (|27|l [as well as the noise correlators in Eqs. 1 -?•!!) 
and l|24|)] only depend on the transmission and reflection 
probabilities T and R. All information about the scat- 
tering phases, from the beam-splitter as well as from the 
propagation in the leads, drops out. As a consequence, 
the correlators are insensitive to dephasing of the orbital 
part of the wavefunction, i.e. processes that cause slow 
and energy independent fluctuations of the scattering 
phases. This insensitivity, different from schemes based 
on orbital entanglement ^2LSi21 can be understood by 
considering the first quantized version [in Eq. (|llfl ] of 
the wavefunction \^i). Any orbital phase picked up by 
an electron in e.g. lead 1 just gives rise to an overall 
phasefactor of the total orbital wavefunction, since each 
term in the wavefunction corresponds to one electron in 
lead 1 and one in lead 2. Moreover, any orbital "pseudo 
spin flip" would imply a scattering of particles between 
the leads 1 and 2 and is not allowed in the non-local 
geometry. 

The situation is different for spin decoherence, energy 
independent spin-flip or spin-dephasing processes tend- 
ing to randomize the spin directions. Spin decoherence 
generally modifies the Fano factors in Eqs. (|27|) and 
(|28|l . Formally, the (mixed) state in the presence of de- 
coherence is described by a density matrix p. Writing 
p in a spin singlet-triplet basis, as shown in the Ap- 
pendix, only the diagonal elements pss (singlet) and 
Pt Tq j Pt + t + j Pt_t„ (triplets) contribute to the current 
correlators. As discussed above, all the three spin triplets 
give rise to the same Fano factors. The spin-triplet 
Fano factors are given by the spin-singlet ones in Eqs. 
(|27|l and (|28ll by changing the sign of the last term 
2RT\H(e\ — £2)| 2 , i-e. from bunching to anti-bunching. 
Using that the sum of the diagonal elements of the den- 
sity matrix is one, i.e. p S s + Pt t + Pt+t + + PT-T- = F 
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the effect of spin decoherence is to renormalize only the 
part of the Fano factors dependent on the dot-level ener- 
gies as 

\H(ei - e 2 )\ 2 -> (2 PSS - l)|£T( El - e 2 )\ 2 . (30) 

The renormalization factor is thus the singlet weight mi- 
nus the total triplet weight, pss — (pt t + Pt + t + + 
Pt_t_) = 2/9,55 — 1- This clearly displays how deco- 
herence, reducing the singlet weight and consequently 
increasing the triplet weight, leads to a cross-over at 
pss = 1/2 from a bunching to an anti-bunching behav- 
ior of the noise correlators. For a completely dephascd 
spin state, with an equal mixture of singlets and triplets 
{pss = PT a T a = Pt + t + = pT-T- = 1/4), the renormaliza- 
tion factor 2pss — 1 saturates at the value —1/2. 

We point out that this discussion might be modified 
when considering other types of effects causing decoher- 
ence, such as e.g. inelastic scattering. A more detailed 
investigation (see e.g. Refs. HH), going beyond the scope 
of the paper, is needed to address these issues. 

C. Spin entanglement bound 

In the absence of spin decoherence the spin state of the 
emitted pair is a singlet, a maximally entangled state. 
For finite spin decoherence, this is no longer the case and 
the question arises how to obtain quantitative informa- 
tion about the spin entanglement from the measurements 
of the current correlators. 

We stress that our interest here is the spin entangle- 
ment of p only. However, p contains information about 
the spin part of the state as well as the energy-dependent 
orbital part, the wave-packet structure of the emitted 
pair of electrons. To quantify the spin entanglement, 
one thus has to consider a measurement sensitive to the 
spin part of p only (see Appendix). One such impor- 
tant example is the cross correlators between the cur- 
rents in the leads 1 and 2 (i.e. without beamsplitters). 
It was shown in a related system in Ref. 15] that these 
cross correlators are simply proportional to the proba- 
bility to jointly detect one particle in lead 1 and one 
in 2. The wave-packet property of the emitted pair re- 
sults only in an overall constant multiplying the proba- 
bilities. As a consequence, a Bell Inequality, derived in 
terms of the joint detection probabilities, could be for- 
mulated in terms of zero-frequency cross correlators. In 
the same way, for the superconductor-dot entangler con- 
sidered here, the spin entanglement of the two emitted 
electrons can in principle be tested via a Bell Inequal- 
ity formulated in terms spin current correlators^* The 
situation is different for the beam-splitter setup, where 
the Fano factors in Eqs. H27|) and l|28|) in general de- 
pend on the wave-packet structure via the dot-level de- 
pendent factor \H(si — e 2 ) | 2 , quantifying the overlap of 
the two electrons when colliding. However, for £\ = s 2 
(i.e. maximal overlap, |i?(0)| 2 = 1) the Fano factors are 
independent of the wave-packet structure of the emitted 



electrons and thus only sensitive to the spin part of p [see 
Eq. EU]. 

The spin part of p can be described by the 4x4 spin 
density matrix p a , rigorously defined in the Appendix 
(note that p, due to the continous energy variable, is in- 
finite dimensional). Formally, p a is the density matrix 
obtained when tracing p, for aligned dot-levels £\ = £ 2 , 
over energies. The question is thus how to determine the 
entanglement of p a . In general, knowledge of all the ma- 
trix elements is needed. This information can however 
not be obtained within our approach, since the Fano fac- 
tors only provides information of the spin singlet weight, 
as is clear from Eq. (|30() . It is nevertheless possible, as 
described in detail in the Appendix, to follow the ideas 
of Burkard and Losses and obtain a lower bound for the 
spin entanglement. 

There are several different measures of entanglement 
for the mixed state of two coupled spin-1/2 systems. Here 
we consider the concurrence*** C, with C = (C = 1) 
for an unentangled (maximally entangled) state. To es- 
tablish the lower bound, it can be shown that the con- 
currence C(p a ) is always larger than or equal to the con- 
currence C(pw) of the so Werner stated described by 
the density matrix pw ■ The Werner state, defined as 
the average of p a over identical and local random ro- 
tations, has the same singlet weight pss as p„. The 
concurrence of the Werner state has the appealing prop- 
erty that it is a function of the spin singlet weight only, 
C w = max{2p ss - 1, 0}. 

The findings above thus lead to the simple and im- 
portant result that the renormalization, Eq. (|3U[1 . of the 
Fano factors in Eqs. (|27|l and 128(1 due to spin decoher- 
ence can be written as (for Cw > 0) 

\H( £l ~e 2 )\ 2 ^C w \H(e 1 -e 2 )\ 2 (31) 

where Cw thus provides a lower bound for the spin en- 
tanglement of the emitted pair of electrons (for the pure 
singlet pss — lj Cw and C(p tT ) are equal and maximal). 
Thus, as long as the Fano factors display a bunching be- 
havior, the spin entanglement is finite, Cw > 0. For a 
cross-over to anti-bunching behavior, Cw — and one 
can no longer conclude anything about the entanglement 
of the spins state. The value of Cw can be extracted 
directly from the experimentally determined Fano fac- 
tors, as the amplitude of the modulation of the Fano 
factors with respect to dot level amplitudes £\ — e 2 di- 
vided by 2RT. The values of R, T can be extracted inde- 
pendently from the Fano factors at dot levels such that 
H( £l - e 2 ) « 0. 

The result in Eq. I|31|) thus provides a simple relation 
between the Fano-factors and the minimum spin entan- 
glement Cw- It is clear, however, that since the Fano 
factors only provide information of the singlet weight, 
full information of the spin entanglement can not be ob- 
tained by the beam-splitter approach employed here. It 
should be noted that the result in Eq. I|31|l is quanti- 
tatively different from what was obtained in Ref. |39| . 
a consequence of the different states considered for the 
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emitted electrons, as discussed above in connection with 
the current correlators. 



VI. TUNNELING THROUGH THE SAME DOT. 

We then turn to the situation when the two electrons 
tunnel through the same dot. To be able to distinguish 
this process II from process I, it is important to study 
the current as well as the noise in detail. The averaged 
current in Eq. 1|19[) . evaluated with the state in Eq. 1)150. 
becomes for reservoirs A and B 

#=f J^dEiRlB^Etf+TlBzm 2 ), 

I'b' = x r dE (m(£)| 2 + ^ 2 (£)| 2 ) . (32) 

It J-eV 

Since the two resonances £\ and e 2 are well within the 
voltage range, i.e. eV — \sx\, eV — |<S2 1 ^> 7, we can 
perform the integrals and get the current^ 

I" = ^7r|S | 2 / 7 , (33) 

the same for both reservoirs a = A,B. We note that the 
two-particle resonance in the current, present in the pair- 
splitting case I, is absent due to the Coulomb blockade, 
as pointed out in Ref . [3^J . A difference from Ref . j32 is 
however that due to the absence of back-scattering at the 
beam-splitter, there is no scattering-phase dependence 
of the current. Consequently, there is no dependence 
on a possible difference in the superconducting phase at 
the two emission points or an Aharonov-Bohm phase 61 
due to a magnetic flux in the area between the super- 
conductor, the dots and the beam-splitter. It should be 
pointed out that this is not a generic result for normal- 
superconducting systems. In a situation with backscat- 
tering, which is inevitable in e.g. the three-terminal fork- 
like geometries, Andreev interferometers, studied exten- 
sively in both diffusive^ and ballistic 63 conductors, the 
current is indeed sensitive to a superconducting phase 
difference as well as the Aharonov-Bohm phase. 

Regarding the spin entanglement, just as for process 
I, no information is provided by the averaged current. 
The same result would have been obtained considering 
an incoherent superposition of two electrons in lead 1 
and two in lead 2, independent on spin state. Turning to 
the current correlators, inserting the expression for the 
state \^>n) into Eq. l|20)) . one gets the expressions for the 
auto-correlations 

Ap 2 f 

S aa = — J dEiRil + R^E)] 2 

+T(1 + T)\B 2 {E)\ 2 + 2Rc [(r*t r ) 2 B*(E)B 2 (E)] } , 

S^ = y / dE{T(l + T)\B 1 (E)\ 2 

+R(1 + R)\B 2 (E)\ 2 + 2Re [(r*t') 2 B* 1 (E)B 2 (E)] } 

(34) 



with Rc[..] denoting the real part, as well as the cross- 
correlations 

S& = S& = x/ dE { RT (\Bi(E)\ 2 + \B 2 (E)\ 2 ) 
-2Re [(r*t') 2 B* {E)B 2 (E)] } . (35) 

The integrals over \Bj(E)\ 2 were carried out above, Eq. 
(|33|l . Performing the integral over B\(E)B2(—E) in the 
limit eV — \ei\,eV — \e 2 \ ^> 7, we get 

/ dEB* (E)B 2 (E) = m \ B °\* . (36) 
J ei - £2 + H 

The expressions for the correlators above give that the 
total noise S 11 of the current flowing out of the super- 
conductors is, 

S" = + SMa + + S'Jb = MIa + Ib), (37) 

twice the Poissonian, describing, just as in case I, an 
uncorrelated emission of pairs of electrons. 

We note that in contrast to the current and the trans- 
port properties in case I, when the two electrons tun- 
nel through different dots, the noise contains information 
about the scattering phases (via r*t'). Quite generally, 
one can write 

(r*t') 2 = RTe^, (38) 

where (f> is a scattering phase of the beam-splitter. Scat- 
tering phases picked up during propagation in the leads 
simply add to 4>. As a consequence, (f> can be modu- 
lated by e.g. an electrostatic gate changing the length of 
the lead 1 or 2 or by an Aharonov-Bohm flux threading 
the region between the dots, the superconductor and the 
beam-splitter. An important consequence of this phase 
dependence of the current correlators is that it can be 
used to distinguish between tunneling via process II and 
between process I, since the current correlators of the lat- 
ter show no phase dependence. This was pointed out in 
Ref. H2. 

This phase dependence shows that the correlators in 
Eqs. (|34|l and 1|35|) are sensitive to dephasing affecting the 
orbital part of the wavefunction. For complete dephas- 
ing, the last term in Eqs. (p^ffl and PBf are suppressed. 
The orbital entanglement in Eq. (|15f) . the linear superpo- 
sition of states corresponding to tunneling through dot 
1 and 2, is lost. This sensitivity to orbital dephasing 
is different from the one for process I discussed above. 
However, again in contrast to process I, the current cor- 
relators are insensitive to spin-dephasing. This can be 
understood by considering the first quantized wavefunc- 
tion l^u), discussed below Eq. (|TS|l . keeping in mind 
that the wavefunction is a direct product of a spin part 
and an orbital part. The spin wavefunction is a singlet, 
xjiX'l ~ x\iX]ji but the orbital wavefunction is a combi- 
nation of triplets, A 1 A 1 + X^X 2 for e\ — e 2 . Since no 
scattering between the leads is possible, i.e. no "pseudo 
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spin flip" , orbital dephasing can not change the triplet 
character of the orbital wavefunction and as a result, the 
spin wavefunction is bound to be a singlet. Thus, the spin 
entanglement in \^fu) is protected against decoherence. 

Turning to the Fano factor gives, the auto and cross 
correlations are 



w 11 

F BB 



l+T 2 + R 2 

2RTRe [e^Hfa - e 2 )] 



and 



*AB — *BA 



= 2RT 
- 2.RTRe [e^H(ei - e 2 )] 

respectively, where H{e\ — £2) is given in Eq. (|29[1 . 



(39) 



(40) 



shows a dip for aligned resonant levels. At E\ — £2 = 0, 
the Fano factor is zero, independent on the beamsplitter 
transparency T. This is a sibnature of perfect bunch- 
ing. For finite phase-difference 0^0, the Fano factor 
becomes asymmetric in S\ — £2, showing a Fano-shaped 
resonance, with the minimum shifted away from e\ = £2. 

The Fano factor for the auto-correlations, for cf> = 0, 
show a corresponding peak for aligned resonant levels, 
reaching 2 for e\ = £2- For finite phase-difference <f> =f 0, 
the Fano factor becomes asymmetric, with the maximum 
Fano- factor shifted away from £\ = £2- 

We point out that similarly to case I, the integrand 
of the cross correlators can be understood by consider- 
ing the basic two-particle scattering processes. They are 
shown in Fig. [SJ the general explanation is along the 
same line as for process I, discussed above. 
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FIG. 6: The Fano factor for the cross-correlations Fab = 
Fba for phase difference <f> = (left panel) cj> — n/2 (right 
panel) and as a function of the normalized energy difference 
(£1 — £2)/7 for various beam-splitter transparencies. 




FIG. 8: Elementary scattering processes (shown at the beam- 
splitter) contributing to the cross correlators Sab- The two 
processes (a) and (b) transport a pair of electrons | |, E) and 
I l,—E) from the superconductor to the reservoirs A and B 
respectively. 



VII. DISCUSSION AND CONCLUSIONS. 



The Fano factor as a function of energy difference 
£1 — £2 is plotted in Figs. [7] and for several values of 
the transparency of the beam splitter. For zero phase dif- 




FIG. 7: The Fano factor for the auto-correlations Faa = 
Fbb for phase difference = (left panel) cj> — n/2 (right 
panel) and as a function of the normalized energy difference 
(£1 — £2)/7 for various beam-splitter transparencies. 

ference = 0, the Fano factor for the cross-correlations 



In conclusion, we have investigated the spin entangle- 
ment in the superconductor-quantum dot system pro- 
posed by Recher, Sukhorukov and Loss |32|. Using a 
formal scattering theory we have calculated the wave- 
function of the electrons emitted by the entangler and 
found that it is a superposition of spin-singlets at differ- 
ent energies, a two particle wavepacket. Both the wave- 
function for the two electrons tunneling through different 
dots, creating the desired nonlocal EPR-pair, as well as 
the wavefunction for the two electrons tunneling through 
the same dot, were calculated. 

The two electrons in the emitted pair collide in a beam- 
splitter before exiting into normal reservoirs. Due to 
the symmetrical orbital state, a consequence of the anti- 
symmetrical singlet spin-state, the electrons tunneling 
through different dots show a tendency to bunch. This 
bunching can be detected via the current correlations. It 
was found that the amount of bunching depends on the 
position of the single particle levels in the dots as well as 
on the scattering properties of the beam splitter. Impor- 
tantly the magnitude of the bunching dependent term in 
the cross correlations was found to be of the same order 
as the bunching independent term, implying that an ex- 
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perimental detection of the bunching, and thus indirectly 
the spin-singlet entanglement, is feasible. 

The current correlators for electrons tunneling through 
different dots were found to be insensitive to orbital de- 
phasing. Spin dephasing, on the contrary, tends to ran- 
domize the spin state, leading to a mixed spin-state with 
a finite fraction of triplets. Since singlet and triplet spin 
states give rise to a bunching and anti-bunching behavior 
respectively, when colliding in the beam-splitter, strong 
dephasing will suppress the bunching behavior and even- 
tually cause a crossover to anti-bunching. To quantify 
the entanglement in the presence of spin dephasing, we 
have derived an expression for the concurrence in terms 
of the Fano factors. In addition, via the current correla- 
tions, it is not possible to distinguish between entangled 
and non-entangled spin-triplet states, since all triplets 
show the same bunching behavior. This implies that the 
method of detecting spin entanglement via current corre- 
lations in the beam-splitter geometry has a fundamental 
limitation compared to the experimentally more involved 
Bell Inequality test. 

We have also investigated the current correlations in 
the case where the two electrons tunnel through the same 
dot. The wavefunction was found to be a linear superpo- 
sition of states for the pair tunneling through dots 1 and 
2. The cross- and auto correlators are sensitive to the po- 
sition of the single-particle levels in the dots, however in 
a different way than for tunneling through different dots. 
Moreover, the correlators were found to be dependent 
on the scattering phases, providing a way to distinguish 
between the two tunneling processes by modulating the 
phase. 
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APPENDIX A 

In the presence of spin decoherence, the state of the 
pair of electrons emitted through different dots can be 
described by a density matrix p, which can be written as 



{q} = {S,T ,T + ,T_},le. 



fdE\A(E)A Y,PW J dEdE ' 
9,9' 

x A(E)A*(E')\%(E)){* q ,(E>)\ 



(Al) 



\*t (E)) 

\* T+ ( E )) 
\*tAE)) 



1 



V2 
1 

71 



6^(^)6^ (-£) - 6 t (^) 6 t (.^1 | ) 



bl^E)bl l (-E) + b\ l (E)b^(-E)\ |0) 
b^(E)b^(-E)\0) 

&l ; (£)4 ; (-£)|0> (A2) 

The coefficients p qq i depend in general on the nature and 
the strength of the spin decoherence. As pointed out 
in the text, only energy independent spin decoherence 
is considered, and consequently the coefficients p qq i are 
independent on energy. 

The current operators conserve the individual spins. 
As a consequence, the off-diagonal elements of p do not 
contribute to the noise correlators. As discussed in the 
text, all triplets contribute equally to the correlators. 
Since the singlet and triplet states contribute with oppo- 
site sign to the last term in Eqs. H3 and the effect 
of spin-decoherence on the Fano factors can be incorpo- 
rated by renormalizing \H{si — e 2 )| 2 — > (2p ss — l)\H(ei — 
£2)| 2 , with the renormalization factor expressed in terms 
of pss only (using Pss + Pt t + Pt + t + + Pt-T- = 1), the 
weight of the singlet component in p. 

It is a difficult (and in general not analytically 
tractable) problem to evaluate the entanglement of the 
full density matrix, since p contains information about 
both the energy-dependent orbital part of the state as 
well as the spin-part. In particular, due to the continous 
energy variable, the dimension of p is infinite. Here, we 
are however only interested in the spin entanglement of 
p. To determine the spin entanglement one has to con- 
sider measurement schemes where the observables O are 
sensitive only to the spin part of p. Such observables 
satisfy the property 

dEdE 1 'A(E)A* '(£')(* q (E)\0\V q >(E')) 
= (* g |Oa|* 9 ') J dE\A(E)\ 2 (A3) 

where \^ q ') are given from \^f q i(E)) in Eq. ilA2() by re- 
moving the energy dependence, e.g. \^t + ) — b^b^JQ). 
The operator O a is a function of the energy independent 
^-operators. Using the property in Eq. (|A3|) we can write 

(O) = tr[pO] = (J dE\A(E)A £ Pw , 

9,9' 

x f dEdE'A(E)A*(E , )(V q (E)\0\y q .(E')) 



(A4) 



noting that the normalization gives ^2 q p qq = 1. The 
index q runs over the states in the singlet-triplet basis 



= J2p qq >(^ q \O a \^ q >) = tr[pM. 

9,9' 

The 4x4 spin density matrix p a is thus 

P* =I^P99'I*9>(*9'I- 
</■'/ 



(A5) 
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It is straightforward to show that for the special p for 
aligned dot levels E\ — e 2 , the current correlators in Eq. 
(121 H are insensitive to the wave-packet structure of p. 
In this case, p a is directly obtained from p by tracing 
over energies. More generally, independent of £i,£2, the 
spin current correlators between lead 1 and 2 (i.e. in 
the absence of the beamsplitter) are insensitive to the 
wave-packet structure of p. These latter correlators can 
be use d to test a Bell Inequality, along the lines of Ref. 

|T3ifl| . 

Our interest is thus to investigate the entanglement of 
p a , conveniently expressed in terms of the concurrence^ 
The concurrence C is defined as 

C(p CT ) = max{o,v / Ai"-VA^-v / ^-V / Al} (A6) 

where the A^s are the real and positive eigenvalues, in 
decreasing order, of papa- The matrix p a is defined as 

Pa = (Cy ® Vy)P*a i. a y ® °v) ( A7 ) 

where a v are Pauli matrices, rotating locally the spins in 
lead 1 and 2 respectively. Importantly, in Eq. \A7\ . the 
density matrix p a is written in the spin-up/spin-down 
basis, i.e. b\^b\ JO) etc. The concurrence is C — for 
an unentangled state and C = 1 for a state which is 
maximally entangled. 

To determine C(p a ), full information of p a is needed. 
In the approach taken here, investigating the spin entan- 
glement via a beam-splitter and current correlators, one 
can however not determine all elements of the density 
matrix p a • As a consequence, the spin entanglement of 
the emitted pair can not be determined precisely. It is 
nevertheless possible, following the ideas of Burkard and 
Loss^S to obtain a lower bound for the spin entangle- 
ment. 

To obtain the lower bound, we first note two im- 
portant properties of C(p a ): (i) C{p a ) is invariant un- 
der local rotations^ i.e. C(p a ) = C{p a ) for p a = 



(t/i ® U<i)Pa(U\ ® U{), where U\ and U 2 are unitary 
2x2 matrices acting locally on the spins in lead 1 
and 2 respectively, (ii) C(p a ) is a convex function, 64 
^iPiC(pi) > C(J2iPiPi), i-e- for a density matrix 
Pa = J2iPiPi> with J2iPi = 1) the entanglement of the 
total density density matrix is smaller than or equal to 
the weighted entanglement of the parts (a consequence 
of information being lost when adding density matrices) . 

Consider then the density matrix pyy obtained by av- 
eraging pa with respect to all possible local rotations 
U ® U , i.e. the same rotation in lead 1 and 2. For- 
mally, pw = ((U ®U)pa(U^ ®U^))u is calculated, where 
(..){/ denotes an average with respect to U, uniformly 
distributed in the group of unitary 2x2 matrices. This 
gives the Werner stated . 

Pw = pss\*s){*s\ + 1 ~£ SS 

X (|* Tn )(* To | + |* T+ )(*T + | + |*T_X*T_|) (A8) 

where we note that the singlet component is unaffected 
by the rotation U ® U . Importantly, the entanglement 
of the Werner state is a function of the singlet coefficient 
pss only. Using the two properties (i) and (ii) of the 
entanglement stated above, we can write 

C(pw) = C[((U<E>U) P a(U^ ®tf))u] 

< (C[(U ®U) P a{U^ ®rf)])v 

= (C( P a))u=C(Pa) (A9) 

This shows that the concurrence of the Werner state 
Cw = C(pw) provides a lower bound for the entan- 
glement of the full spin state C(p a )- The concurrence 
of the Werner state is Cw — rnax{2p5s — 1,0}. The 
renormalization of the Fano factors in Eqs. (|27|l and 
(|28|l due to spin decoherence can now simply be written 
\H( £l - e 2 )\ 2 -> C w \H(ei - e 2 )\ 2 where C w > is a 
lower bound for the concurrence of the spin state in the 
presence of decoherence. This is Eq. (|31f) in the text. 
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